Positive cross-correlations induced by ferromagnetic contacts 
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Due to the Fermionic nature of carriers, correlations between electric currents flowing through 
two different contacts attached to a conductor present a negative sign. Possibility for positive 
cross-correlations has been demonstrated in hybrid normal/superconductor structures under certain 
conditions. In this paper we show that positive cross-correlations can be induced, if not already 
present, in such structures by employing ferromagnetic leads with magnetizations aligned anti- 
parallel to each other. We consider three-terminal hybrid structures and calculate the mean-square 
correlations of current fluctuations as a function of the bias voltage at flnite temperature. 
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I. INTRODUCTION 

Due to the discrete nature of charge carriers, the elec- 
tronic current flowing through a conductor is subject to 
time-dependent fluctuations around its mean value. Cor- 
relations of such fluctuations are of great interest since 
they provide more information, with respect to average 
currents, on the physics of transport. Shot noise, defined 
as the mean-square fluctuations of the current flowing 
through a given terminal at zero temperature, has been 
extensively studied in a wide variety of systems (for a re- 
view on the subject see Ref. [0). Here we are interested 
in the correlations of the current fluctuations between 
different contacts (cross-correlations). Biittiker first 
pointed out that while cross-correlations can either be 
positive or negative for Bosons, for Fermions they are 
necessarily negative, both in the equilibrium situation 
and in the transport regime (see also Ref. Q). 

At this point a question arises naturally: what hap- 
pens in the presence of superconductivity, where charge 
current is carried by the superconducting condensate in- 
stead of by Fermionic excitations. In the sub-gap regime 
the presence of the condensate manifests itself in the dou- 
bling of shot noise, as found for normal/superconductor 
(NS) interfaces theoretically in Refs. and experi- 

mentally confirmed in Refs. |0-|l^. Moreover, in Refs. 
p] |l5| it has been shown that multi-terminal hybrid 
structures containing superconducting inclusions or leads 
can, in some cases, exhibit positive cross-correlations. 
In particular, in Ref. |jll| current correlations have been 
studied as a function of the phase difference between 
the superconducting order parameter of two supercon- 
ducting pieces. In the sub-gap regime, positive cross- 
correlations have been found within a range of </> in the 
case of a floating superconducting electrochemical poten- 
tial. In Refs. [|l2|,^ 15 1 a fork-like geometry with a super- 
conducting leg has been considered and the possibility of 



positive cross-correlations has been found for bias volt- 
ages smaller than the superconducting energy gap. In 
Ref. positive cross-correlations have been predicted 
in a NS heterostructure with a weakly coupled tip. It 
has been argued, however, that such positive correlations 
vanish in the limit of large number of channels or in the 
presence of disorder. 

It is worthwhile mentioning that all the above results 
are relative to a fully phase-coherent transport regime. 
A doubling of shot noise, however, has been also found 
within a semi-classical theory p^ , p^ , therefore proving 
that phase-coherence is not required for such an effect to 
occur. More importantly, in hybrid systems in which a 
semi-classical treatment is valid, cross-correlations have 
been demonstrated to be always negative. 

Recently, advances in nano-fabrication techniques of 
hybrid normal- metal/ferroniagnet structures has led to 
the advent of spintronics [|l8|. In such a context the 
spin degree of freedom of electrons is exploited (instead 
of their charge) and spin-current correlations can be ac- 
cessed. In this paper we introduce this new ingredient 
and we address the question whether ferromagnetism has 
an effect on the sign of the cross-correlations. 

As far as transport is concerned, a ferromagnet can be 
thought as a conductor in which spin-up and spin-down 
electrons contribute by a different amount to the total 
current (the spin degeneracy is broken). For example, in 
the limit of a 100% polarized ferromagnet (the so called 
"half metal"), current is carried by a single spin species. 
Due to the fact that in a s-wave superconductor spin-up 
particles are coupled to spin-down holes, non-trivial ef- 
fects are expected on the cross-correlations between two 
ferromagnetic contacts having magnetizations in oppo- 
site directions. We shall show that a positive value of 
cross-terminal correlations in such hybrid systems can be 
induced (as noticed in Ref. |jl^) or, if already present, en- 
hanced. Remarkably, positive cross-correlations are pre- 
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dieted over a wide range of bias voltages, even above the 
energy gap, for structures which exhibit negative corre- 
lations when the ferromagnets are replaced with normal 
metals. As expected, such an effect gets stronger as the 
ferromagnetic exchange field is increased and reaches a 
maximum for exchange fields corresponding to the limit- 
ing case of half-metallic leads. The effect remains essen- 
tially unchanged even in the presence of disorder, due to 
impurities and lattice imperfections. 



order parameter, non-zero only in the superconducting 
regions. Let us consider, for example, a unit flux of par- 
ticles originating from the j-th reservoir. The scattering 
state for the j-th lead (assuming for simplicity a single 
open channel per lead) takes the form 



u{x) 
v{x) 



(2) 



II. THE MODEL 

We consider a realistic hybrid 2D system attached to 
three multi-channel terminals. As sketched in Fig. |l|, it 
consists of a superconducting island (which defines the 
scattering region) connected to one normal metallic lead 
(1) on the left-hand-side and two ferromagnetic leads (2 
and 3) on the right-hand-side (similar structures have al- 
ready been considered in Refs. [|9-22|). The directions 
of magnetization of leads 2 and 3 are chosen to be anti- 
parallel to each other. The two ferromagnetic electrodes 
are kept at the same potential with respect to the normal 
lead {Vi — V2 = Vi — V3 = V). For definiteness we as- 
sume the magnetization of lead 2 directed upward while 
the magnetization of lead 3 directed downward. 

Shot noise and cross-correlations are calculated within 
the multiple scattering approach. The scattering matrix 
S completely characterizes a given structure (scatterer) 
attached, through perfect conductors (leads), to N reser- 
voirs of particles. S is defined via the asymptotic wave 
functions of the leads, known as scattering states. When 
the scatterer contains superconducting regions, the scat- 
tering states are solutions of the Bogoliubov-de Gennes 
equation |23 



Hp A 
A* Hh 



u{x) 
v{x) 



= E 



u{x) 
v{x) 



(1) 



In Eq.(^ Hp is the Hamiltonian relative to the parti- 
cle degrees of freedom, Hh is the Hamiltonian relative 
to hole degrees of freedom and A is the superconducting 



whereas for the i-th lead, with i ^ j , it reads 



u{x) 
v{x) 



(to) 



(3) 



Here (t'o)^ and {ra)j are, respectively, normal and An- 
dreev reflection amplitudes, while (to)ij a-nd {ta)ij are, re- 
spectively, normal and Andreev transmission amplitudes. 
In Eqs. (||) and (^) Vj {ujj) is the group velocity in lead 
j for particles (holes) and kj {qj) is the wave- vector for 
particles (holes) . S is the matrix of scattering amplitudes 
of dimension {2N x 2A^) defined by 



I out) = S I in) 



(4) 



which relates the amplitudes of an incoming wave func- 
tion ( I in) ) to the amplitudes of an outgoing wave func- 
tion ( I out) ) . We assume the convention that the upper 
N elements of such vectors are relative to particles and 
the lower ones to holes. 

Following Refs. [p]J^j2^ the noise power su is defined 
as the Fourier transform of the mean correlations of the 
current fluctuations in leads i and I: 
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{Ah{t)AIi{t') + AIi{t')Ah{t)) 



(5) 



where AIi{t) = ii{t) — (li) is the current fluctuation op- 
erator in lead i at time t and the angle brackets (. . .) 
denotes a statistical average. By defining the time- 
dependent current operator as 
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(6) 



3/31) 
j'0'b- 



one obtains the following expression for the zero-frequency noise power 

=2 /-oo 



2g2 poo ^ 

su{V,T) = — dEj2Yl A'^p,.^^,p,,,ii,E,E)A}p,,,.^^p,{l,E,E) {f,p{E) [1 - frp>{E)\ + J,.p.{E) [1 - J,p{E)\} 

aa' jfb 



(7) 

In Eq. (|^) o}jpij{E) is the creation operator for an incoming particle (/3 = +1) or hole {f3 ~ —1) with open channel 
index b at energy E originating from lead j. The quantity A°:pj^.-,p,^,{i,E,E') is defined in terms of the scattering 
matrix as 
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(8) 



whereas fjp^E) is the Fermi distribution function for 
particles {(3 = +1) or holes {f3 = —1) at energy E in 
reservoir j: fj(3{E) = [1 + exp(£; - (3{eVj - fi)/kT)]-\ 
H being the condensate electrochemical potential. Note 
that S^iaa),{ji3b) is the scattering amplitude for a particle 
(/? = +1) or a hole {f3 = —1) originating in lead j and 
relative to channel b to be transferred into lead i, channel 
a as a particle (a — +1) or a hole {a = —1). 

To evaluate the elements of the S-matrix we make use 
of an efficient recursive Green's function technique ||2^ ] 
which is based on the tight-binding model. We there- 
fore describe the system under investigation in terms of 
a 2-dimensional tight-binding Hamiltonian on a square 
lattice: 



r,k,a,l3 



(9) 



where is either a particle (a = +1) or a hole 

(a = — 1) degree of freedom of an atomic-like state (Wan- 
nier function) centered at a site I. In the case of a perfect 
crystal without imperfections we set the on-site energies 
as follow 



a(eo - h) 



(10) 



where ep is a real constant determining the position of 
the Fermi energy and h, according to the Stoner model, 
is the ferromagnetic exchange field, non-zero only in leads 
2 and 3 (positive in lead 2 and negative in lead 3). Note 
that the spin degrees of freedom are not coupled, since 
the magnetizations in leads 2 and 3 are assumed to lie 
along the same direction. For the hopping potentials we 
have 



with I and k first nearest neighbours, and 



V* 
i,- 



.1,+ 



(11) 



(12) 



since in local s-wave superconductors the particle-hole 
coupling is on-site. Note that while F is a real constant 
determining the band width, the superconducting order 
parameter A is complex. The presence of impurities and 
lattice imperfections is modeled as Anderson disorder, 
whereby eo in (^0|) is replaced by a uniformly distributed 
random number in the range [eo — W/2,eo + W/2], W 
being the disorder strength. 

The total Hamiltonian is block diagonal so that one 
can separately consider the Hamiltonian relative to spin- 
up particles and spin-down holes from the Hamiltonian 



relative to spin-down particles and spin-up holes. Since 
they are decoupled, no correlations are present between 
spin-up and spin-down particles (or holes) . From this fol- 
lows that, as long as the system is symmetrical under the 
exchange of leads 2 and 3, the results for spin-down po- 
larized injected current are equal to the results for spin- 
up polarized injected current once indices 2 and 3 are 
exchanged. For simplicity, then, we assume the current 
injected from lead 1 to be 100% spin-up polarized and for 
this reason we can drop the spin index in the definition 
of the on-site energy. In the recursive Green's function 
technique the scattering amplitudes are calculated nu- 
merically from the total Green's function of the system, 
which in turn is evaluated through the effective Hamil- 
tonian of the scattering region and the Green's functions 
of the isolated leads iGSl . 



III. RESULTS 

In our calculations the scatterer (shaded region in Fig. 
|l|) is a superconductor 15 sites long (in unit of lattice 
spacings) and 15 sites wide. The leads are 3 sites wide: 
lead 1 is attached in the central position on the left- 
hand-side, whereas leads 2 and 3 are placed symmetri- 
cally with respect to the center of the right-hand-side at 
a distance of 5. In arbitrary units, the order parameter is 
set to A = 0.1 uniformly in the superconducting region, 
therefore neglecting the suppression of the superconduct- 
ing order parameter in the regions adjacent to the ferro- 
magnets. Such a suppression would lead to very small 
corrections, since it is expected to occur over very short 
lengths (units of nanometers for hard ferromagnets) , or- 
ders of magnitude smaller than the superconducting co- 
herence length. The temperature in the reservoirs is set 
to T = 0.01, while the Fermi energy is about 1.6. With 
these values we have three open channels in the normal 
lead for particles and holes at the Fermi energy and a su- 
perconducting coherence length of the order of 10 lattice 
constants. Note furthermore that the condensate elec- 
trochemical potential is set to /x = 0. 

As a reference situation we first consider the case where 
leads 2 and 3 are normal-metallic {h = 0). In Fig. ^ and 
Fig. H we plot the noise power relative to the different 
pairs of leads as a function of the bias voltage V. When 
the superconductor is in the normal state (Fig. |2|) we 
find, as expected, that the shot noise in lead 1 is posi- 
tive whereas the cross-correlations are negative. In the 
presence of superconductivity (Fig. ^), even though the 
magnitude of the cross-correlations changes, the sign re- 
mains negative. We now arrive at the main result of the 
paper, namely, that a positive sign can be induced by re- 
placing leads 2 and 3 with ferromagnets characterized by 
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magnetizations in opposite directions. This is shown in 
Fig. ^ for the limiting case of half-metals, where cross- 
correlations exhibit a positive sign. While si2 changes 
sign as V is increased, the effect on S23 is more dramatic 
as it maintains a positive sign over the whole voltage 
range, with a large magnitude. It is remarkable that 
S23 remains positive over wide ranges of voltages, even 
exceeding the superconducting energy gap. Of course 
S12 and Si3 are less affected by ferromagnetism than S23, 
since lead 1 is normal-metallic. In the following we will 
concentrate on the behaviour of S23- 

In Fig. ^ we plot S23 as a function of the magnitude 
of the exchange field h in the ferromagnetic leads 2 and 
3 for four different values of bias voltages. First note 
that the four curves show the same behaviour, more pro- 
nounced in magnitude for larger V. For zero exchange 
field, i.e. with normal contacts, S23 takes negative val- 
ues despite the presence of superconductivity. For finite 
values of h, S23 increases turning to positive values, and 
thereafter increasing further rapidly reaching a plateau. 
The plateau corresponds to values of exchange field rela- 
tive to full polarization, therefore it is an artifact of the 
small number of open channels. 

As already noted in Refs. the reason for which 

positive cross-correlations arise in the presence of super- 
conductivity can be understood by expressing the current 
operator of Eq. as a sum of particle and hole contri- 
butions: 

m = i+{t)+ir{t) . (13) 

By substituting Eq. ( p^ in Eq. (||), the noise power su 
[Eq. (^] can be written as the following sum: 

s^l - 4+ + sr + 4" + sTi^ (14) 

where sj"*" is the particle-particle (pp) cross-correlation, 
sfi~ is the particle-hole (ph) cross-correlation and so on. 
Quite generally it was shown in Rcf. that while 4^ 
and s~i~ are negative, 4~ ^-^^ ^ii^ take a positive sign. 
Intuitively, one can understand that pp (and hh) cor- 
relations are negative because of the anti-bunching be- 
haviour due to the Fermionic nature of electrons 
For the same reason, ph and hp correlations present a 
positive sign simply because particle and hole have oppo- 
site charge (note that when the superconductor is in the 
normal state ph and hp correlations are zero). The sign 
of the overall cross-correlation is therefore determined by 
the competition between the positive (ph and hp) and 
the negative (pp and hh) contributions. The important 
observation is that pp and hh cross-correlations can be 
drastically suppressed in the presence of spin-polarized 
transport, for example by allowing only spin- up polar- 
ized current in one lead and only spin-down polarized 
current in the other. As a result an overall positive sign 
will be induced. 

In the system under investigation the different contri- 
butions to S23, in the zero-temperature, small bias volt- 
age limit and with fi set to zero, take the simple form 



= -("7)-^ [T2^^T^^'■ + T^rT^-] , (15) 

where T°l^ is the transmission probability, evaluated at 
E = 0, for a particle (/3 = -1-1) or a hole (/3 = — 1) 
injected from lead 1 to be transmitted into lead j as a 
particle (a = -1-1) or a hole (a — —1). In the limiting 
situation where leads 2 and 3 are half-metallic with mag- 
netizations aligned anti-parallel to each other (for exam- 
ple upward for 2 and downward for 3), we have that only 
positive correlations are non-zero. Particle-particle and 
hole-hole correlations, in fact, vanish (3^3^ = 5^3" = 0) 
because transmission of spin-up particles is forbidden in 
lead 3 and transmission of spin-down particles is forbid- 
den in lead 2: Tgtj" = T^i" = 0. For weaker ferromag- 
nets 523* finite, but nevertheless suppressed with re- 
spect to ph correlations so that S23 remains positive. It 
is worthwhile noting that it is possible to switch the sign 
of S23 from positive to negative by changing the relative 
alignment of the magnetizations in the two ferromagnetic 
leads from anti-parallel to parallel. In the latter case, in 
fact, Tg^" and T2^" in Eq. ( |l5| ) are both equal to zero 
in the limit of half- metallic leads, so that only 5^3^ is 
left finite, analogously to what happens when the super- 
conductor is turned normal. For completeness, note that 
by summing over the indices a and 7 in Eq. ( [Tsl ) and 
employing the particle-hole symmetry one obtains the 
simple expression 

l23 = -^(|(i0)31p-|(ia)31P) X 

x(|(to)21p-|(ta)21P) . (16) 

We can also understand the reason why positive corre- 
lations in Fig. ^ survive over such a wide voltage range. 
While pp correlations are zero, ph correlations tend to 
vanish only for quasi-particles energy well above the su- 
perconducting gap (about one order of magnitude), i.e. 
when Andreev processes are negligible. Furthermore, the 
effect on the sign of S12 and S13 is less pronounced than 
on S23 , because in the former cases pp or hh correlations 
can still be large. It is worthwhile mentioning that the 
presence of barriers at the interfaces between electrodes 
and the superconductor is expected to affect the cross- 
correlations (analogously to the findings of Ref. p3[). 
Eq. ( p^ shows that the sign of S23 can change with 
the barrier strength, since such a sign is determined by 
the interplay between normal and Andreev transmissions 
between leads 1 and 2, and 1 and 3. In the case of half- 
metallic contacts, it is clear that cross-correlations S23 
will still be positive regardless of the value of the barrier 
strength, because negative contributions (pp and hh cor- 
relations) are always zero. However, we expect that the 
magnitude of S23 might depend on the barrier strength 
non-monotonically. In fact, while normal transmission 
coefficients decrease with increasing barrier strength, An- 
dreev transmissions do not monotonically depend on the 
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barrier strength, but instead first increase and then de- 
crease. 

To conclude, we have also checked that the presence of 
weak disorder within the superconductor, due to impuri- 
ties or lattice imperfections, does not significantly change 
the results. Using the Anderson model described above, 
we have considered the same system studied in the clean 
limit for different values of the disorder strength W. Even 
without performing a systematic study of disordered sys- 
tems, we have found that S23 always remains positive over 
the whole voltage range, though the shape of the curve 
changes. On the contrary, for some disorder realizations 
we have found that S12 and S13 do not change sign re- 
maining negative over the voltage range considered. For 
the future it would be of interest a systematical study of 
disordered systems in the different transport regimes. We 
argue that even in the diffusive regime and with a large 
number of open channels in the contacts such positive 
sign survives. 



IV. SUMMARY 



The study of electronic current-current correlations in 
mesoscopic structures has been receiving interest since 
the last decade. On the one hand, it has been found 
that the correlations between different leads are always 
negative. On the other, in hybrid NS structures it has 
been proved that such correlations can, in some cases, 
be positive. Here we have considered the sign of cross- 
correlations in a hybrid NS three-terminal structure in 
which ferromagnetic leads are employed. We have found 
that a positive sign of the cross-correlations between two 
ferromagnetic leads is induced when their magnetizations 
are aligned anti-parallel to each other. It is remarkable 
that such a positive sign persists over a wide range of 
voltages exceeding the value corresponding to the su- 
perconducting energy gap. Furthermore positive cross- 
correlations have also been found between normal metal- 
lic and ferromagnetic leads. The origin of this effect has 
been attributed to the suppression of particle-particle 
and hole-hole cross-correlations (which contribute with 
a negative sign) in favour of particle-hole correlations 
(which contribute with a positive sign), due to the strong 
spin-asymmetry introduced by ferromagnetism. We have 
moreover realized that it is possible to switch the sign of 
S23 from positive to negative by changing the the relative 
alignment of the magnetizations in the two ferromagnetic 
leads from anti-parallel to parallel. Such effects are ro- 
bust against the presence of weak disorder within the 
superconducting region. To conclude, we would like to 
enphasize that our approach can be applied to disordered 
diffusive samples as well as clean and ballistic ones and 
to any geometry. Furthermore, it can be easily extended 
to finite frequencies. 
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FIG. 1. The system consists of a superconductor (S) con- 
nected to 3 contacts. Terminals 2 and 3 are made of ferro- 
magnets, with anti-parallel magnetizations, while terminal 1 
is normal-metallic. Lead 1 is kept at positive potential Vi with 
respect to leads 2 and 3. The relevant geometric lengths are 
of the same order of the superconducting coherence length. 
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FIG. 2. Noise power as a function of bias voltage V with S 
in the normal state and /i = in leads 2 and 3. Due to the 
geometrical symmetry sis is exactly equal to s\2- 
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FIG. 3. Noise power as a function of bias voltage V with S 
in the superconducting state and /i = in leads 2 and 3. The 
vertical line denotes the voltage corresponding to the super- 
conducting energy gap. Due to the geometrical symmetry S13 
is exactly equal to S12 . 
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FIG. 4. Noise power as a function of bias voltage V with 
S in the superconducting state and half-metallic leads 2 and 
3. Note that cross-correlations are plotted on a different scale 
with respect to s\\. The vertical line denotes the voltage 
corresponding to the superconducting energy gap. 
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FIG. 5. Cross-correlation S23 as a function of the magni- 
tude of the exchange field h in leads 2 and 3 for four different 
bias voltages. 
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